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SYNOPSIS 


GREEN'S FUNCTION METHOD FOR ELECTRON DIFFRACTION 
AND BAND STRUCTURE PROBLH^S 
Vibhakar Vasantrao Bhokare 
Ph.D. 

Department of Physics 
Indian Institute of Technology, Kanpur 
December 1973 


The Green's function method of electronic band structure calculation 
due to Korringa-Kohn-Rostoker is a powerful computational tool which has 
the virtue of rapid convergence of eigenvalues. But it has not found as 
much application as it deserves because a straightforward application of 
the method requires enormous computational efforts in many cases. Methods 
similar to the Green's function method have also been formulated in the 
past for the LEED Qjow-Energy Electron-Diffractio^ problem excluding 
surface effects. But there was no actual numerical computation with such 
methods due to the computational difficulties. The present thesis, 
"Green's Function Method for Electron Diffraction and Band Structure 
Problems", attempts to get actual numerical results both for the electron 
scattering problem and the band structure problem of hep ||\exagonal close- 
packecQ lattices, after proper formulation and careful simplification of 
the Green's function method. The results show that in spite of the 
approximations involved in the variational procedure used, one gets 
reasonably good numerical estimates for the electron scattering problem. 

The results for electronic band structure of Beryllium and Titanium are 
close to those obtained by the APW j^gmented Plane Wav^ method. 
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The problem of elastic scattering of electrons from semi- inf inite 
perfect crystals is considered first. This problem arose from and is 
intimately connected with what is now the vast and complicated subject of 
LEED problem. An approximate variational approach is described for 
cODiputing elastic scattering of electrons from the perfect semi-infinite 
crystals. It is shown that the Green's function method of electronic 
band structure calculation can be conveniently used for such computations 
with the assumption of muffin-tin potential at the lattice sites. 

Numerical computation of the intensity in the (0,0) direction for 
Aluminium crystal is carried out. In this connection, the high-energy 
band structure of Aluminium is also computed. Results for (0,0) intensity 
is compared with some theoretical and experimental results of LEED for 
Aluminium. The variation method described here immensely simplifies the 
calculation and yields fairly good results. It may therefore be used as 
a quick and inexpensive method for computing elastic LEED intensities. 

The problem of computing electronic band structure of hep metals 
by Green's function method is considered next. Although this method has 
been used extensively for the cubic metals, it has found very limited 
application to other structures, mainly due to the computational diffi- 
culties for complex lattices and inadequacy of muffin-tin approximation. 
However, the case of hep metals is Interesting because muffin-tin approxi- 
mation in this case is comparable to that for face-centred cubic metals. 
Therefore, it appeared worth-while to study hep metals by Green's function 
method and exploit its virtue of convergence of eigenvalues with small 
size determinants. The Kohn-Rostoker method for complex crystals is 
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employed and group-theoretical factorization of the secular determinant 
is performed by using the symmetry-adapted hexagonal harmonics, which 
greatly simplified the computation. Using this method the band structures 
of Beryllium and Titanium are computed. With enough experimental data on 
band structure, Beryllium is an interesting test case. Titanium is a 
typical transition metal and with little experimental data and limited 
theoretical calculations, it is expected to exhibit the general features 
of hep metals* The energy-bands, Fermi energies, density-of -states curves, 
and coefficients of electronic specific-heat are computed for both, 
Beryllium and Titanium. The results are compared with the earlier theore- 
tical and experimental results, and they are found to be in good agreement. 
The band structure of Titanium with respect to its exchange parameter, ^ , 
is studied at preliminary level. Although sufficient experimental data 
on its band strutture and Fermi surface are required to definitely fix the 
value of the exchange parameter,' It is concluded that 2/3 j^#C^3/4 Is 
the most interesting range. 

The final conclusion of the study presented in the thesis is that 
with proper formulation and careful simplification, it is possible to 
employ the Green's function method of Korrlnga-Kohn-R-ostoker both for the 
electron scattering problem and the problem of electronic structure of hep 
metals leading to quite satisfactory numerical results. 



CHAPTER I 


INTRODUCTION 

The Green's function method of electronic band structure calcula- 
tion for perfect solids was first proposed by Korringa,^ in 1947. His 
approach was based upon the multiple scattering of electron waves in 
solids l^olid being considered as assembly of spherical scatterer£| . 

Later, the same method was proposed on different lines by Kohn and 
2 

Rostoker, in 1954. A variational procedure was Introduced by them to 
solve the integral form of the Schrodinger ' s equation for an electron 
in a periodic potential. Though these two procedures were independently 
given, they are essentially equivalent and are often called the KKR 
I^Korringa-Kohn-RostokerJ method. 

From a mathematical point of view, it has the virtue of rapid 
convergence of eigenvalues. This is intimately connected with the fact 
that the method deals with the scattering of slow electrons in a periodic 
potential, by employing the expansion of the Bloch waves in spherical 
harmonics. It has proved to be rapidly convergent owing to the well-known 
fact, that the partial wave analysis for scattering by single atoms is 
rapidly convergent for slow electrons. 

Computation-wise the Green's function method is a powerful tool. 

But the method has not fcund as much application as it deserves, mainly 
because a straightf oiward application of the method requires considerable 
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computational efforts. In many cases, proper formulation and careful 
simplification is essential to extract numerical results . Such an 
attempt will be made in the present thesis, whore numerical results 
are obtained for both the electron scattering problem and the band 
structure problem of hep ^hexagonal c lose-packec^ structures. 

Ihe KKR method seems to be very useful for theoretical investiga- 
tion of the low energy electron scattering from crystals because basically 
it is a scattering theoxry. In fact, in the recent years, the KKR type 

fomalisms for LEED Qjow-Energy Electron-Dif f ractioiTI Intensity calcula- 

^7 8 

tions have been given by many workers’^" including the author. The fact 
that intimate connection exists between the electron-diffraction theory 
and the band theory, is very helpful for LEED intensity calculation in 
view of the detailed computational facilities existing for the band 
structure calculation. The study of elastic scattering of the low 
energy electrons from the perfect semi-infinite lattices constitutes 
the first part of this dissertation, where it will be shown that the 
well-known KKR method can be used conveniently for the computation of 
electron-diffraction intensities with the additional assumption of 
muffin-tin potential at the lattice sites. 

Although the KKR method was formulated in 1954 it was put to 

9 

widespread practical use in 1961 by Ham and Segal 1 who modified the 
KKR procedure in several Important respects. After that in the last 
decade, the method has been extensively applied for the energy band 
calculation of the cubic metals. But the method has found limited^^ 
application to structures other than the cubic metals, due to the 
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inadequacy of muffin-tin approximation and the computational difficulties 
for the ccanplex crystals. However, the case of the hep metals is inter- 
esting because muffin-tin approximation in this case is comparable to 
that for the fee |jace-centred cubi^ metals. Therefore, it appears 
worth-while to apply the method to hep metals and exploit its virtue of 
convergence of eigenvalues with small size secular determinants. The 
electronic band structure of hep metals v/ill be described in the later 

12 f- 

part of this work. The KKR formalism for complex crystals Ji.e., 

crystals with more than one atom per unit cel^will be employed and 

13 

group-theoretical factorization of the determinant will be described 
which immensely simplifies the computation. 



CHAPTER II 


THE GREEN ' S FUNCTION METHOD 

Since the comprehensive reviews of the Green's function method 

2 9 10 

are available in the literature, ’ * only a brief account of the 
method is presented in this chapter for the sake of completeness. In 
the following sections, we give definition of Green's functions followed 
by the description of the method. 

2.1 Green's Functions •• Definition : 

Consider the Inhomogeneous equation 

L^(r) = C^(r), ' (2.1) 

/I 

with L a Hermitian differential operator, with 
function of the variable _r» and A a given constant. The function 
'^’^r), which satisfies (2.1), obeys the usual type of the homogeneous 
boundary conditions. 

Formally, then, a particular "solution" of (2.1) is conveniently 
constructed in terms of the Green's function G(r,r' ) for the operator 
[li- X^and the given boundary conditions as: 

A|/‘(r) = J G(r,r‘) C§(r') dr' , 

Over Entire 
Space 


continuous 


( 2 . 2 ) 
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an integral equation fory^ir). 

The Green*s function, defined as that .solution of 

for which C^(r) = (^^(r’^r* ) • Thus the Green* s function satisfies the 
equation 

LG(r,r') - ^G(r,r ' ) = c^(r-r'), (2.3) 

with the same boundary conditions as on the function '^(r) . 

The Green's function is of much theoretical importance because it 
transforms a differential equation with suitable boundary conditions to 
an integral equation. And, the statement of a problem in integral form 
is preferred to that in differential form since the integral form is a 
complete statement of the problem in the sense that it contains the 
boundary conditions associated with the problem. Moreover, sometimes 
there is mathematical convenience and ease with integral equations. 

Now, consider the three-dimensional one-electron Schrodinger 
equation in atomic units ; 

- U(r) + E]^(r) 0, (2.4) 

where U(r) and E are potential energy and total energy, respectively. 
This is the basic equation to be solved for the electron diffraction 
and the band structure problems, with different boundary conditions, 
where U(r) is a potential energy of an electron in a tri-periodic 
crystal , 

Treating, U(r^) ''^^Cr ) ='V(r), in (2.4) temporarily as a given 
inhomogeneity even though it contains the unknown function'Y^(r ) , 


it may be written as 
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+ E]A^(r) = VCr) 


(2.5) 


On comparing (2.5) with (2.1), it becomes obvious that the Green's 
function for the differential operator + eJ will satisfy 


[V^ + EjG(r,r') = ^:f(r-r’ ) , (2.6) 

and the particular solution of (2.5) will be given by an Integral equa- 
tion like (2.2). 


2 .2 The Description of the Method : 

For the band structure problem, one seeks solutions of the 
Schrodinger equation (2,4) for a potential U(r) which is periodic in 
lattice translations, so that the Bloch condition 



is satisfied. Here k is the crystal momentum vector, and r is any 
translation vector of the crystal lattice. Equation (2.4) is in atomic 
units where unit of length is Bohr radius and unit of energy is Rydberg, 

Generally, one is interested in calculating the energy values, 

E = E(k), as a function of arbitrary k lying within the first Brillouin 
zone [b.Z.J of the reciprocal lattice space for which solutions of the 
form of (2.7) exist. 

2 

Kohn-Rostoker proceeded to solve this problem from the equivalent 
formulation that '^|/'(r) is also a solution to the integral -equation (2.2), 
whetTs G(r,r ' ) = G(t;-x' ) satisfies (2.6) and ^^£) ^''VXr). 
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Since U(r) is a periodic potential, 

U(r) Svir-r^), 

s 

where V(r-r ) is the potential at r due to the ion at the s lattice 

■“ ""S ”• 

site, Tg * For *'muff in-tin'* model of potential V(r), that is, 

if V(r) is Spherically symmetric within a sphere of radius r^ centred 
around each ion and inscribed in the unit cell ^Here a "simple" crystal 
is considered with a single ion in each unit cell and a centre of 
inversion at the and if V(r) is constant in the region between the 

spheres with which the zero of the energy is made to coincide; then the 
integral equation 


^(r) = J G(r,r ' ) U(r‘ ) ^(r ' ) dr' 

= S 1 G(r,r') VCr'-r^,) ^(r-) dr' 
s' 

is a sum of contribution from the separate cells* Since the muffin-tin 
spheres do not overlap, this is really speaking a sum of the separate 
Integrals, as may be demonstrated by introducing a variable P such 
that r = r + £ in the cell around the lattice sit^i r , Thus 


= S J V(i')Y^(r^.+il) dll. 


over cell 
vol. r* 

Using Bloch condition (2.7), 


ik*(r ,-r ) 

e s s £g‘£gi^ 




V(f')^(l') dll 


(2.8) 
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where, p (£, £' ) = e ® G(_£, £' + r „), is known as the 

s" ® 

"Structural Green's Function", and r „ = r , -r . Henceforth, the 

' '“S “S ‘**S 

variables £ and £* will be denoted by r and r‘, respectively. ^ ) 

is a function of two variables r and r* which refer to the points in the 
same unit cell* 


The integral equation (2.8) has interesting implication that in 
order to solve the band structure problem, it is sufficient to consider 
the solutions of the Schrodinger equation 


- V(r) + E^-^r) = 0, 

within a single unit cell only. 

The Green's function, (t ,r ' ) , may also be written as 


*^(r,r') = - -rp X 


i(k+K ) ‘(r-r' ) 




(2.9) 


where is the unit cell volume and the sum runs over all reciprocal 
lattice vectors - It is evident from (2.9) that^(r,r') is a function 
of k and E, and satisfies the relations: 


(r ' ,r) = ^*(r ,r ' ) , 


and 


( 2 . 10 ) 


(r-tTg,r' ) = e 


ik*r 


^(r,r' ). 


Kohn-Rostoker introduced a variational procedure leading to 


a stationary value for the energy to solve the integral equation (2,8). 
The Integral equation (2.8) may be derived from the variation principle 
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where 


(2ai) 


r 

7\ = j 'Y^(r) V(r ) ^(r 


) dr 


f I u) V<r)^(r,r')V(r’)y(r')di‘di 
t A/ 


I X 


This variational principle has the important property that ^y\_ vanishes in 
first order for all variations from the solution of (2.8) regardless of 
whether or not the variations satisfy the boundary conditions. 




For proper handling of the singularities of , a limiting procedure 


such as 


/v = 

6^0 £ 


where 


A^= I dr V(r)^'^(r) - j 


(2.12) 

(r,r') V(r«) 


r^r^-2€ 


r'4:r^-£ 


l'^r')dr , 


is used. ^ is a positive quantity, and r^^ is the muffin-tin radius. For 
an arbitrary potential, even this variational procedure would be extremely 
complicated to use. However, the method simplifies enormously for the 
muffin-tin potentials j ^Muffin- tin approximation is precisely the major 
limitation of the method^. If the zero of energy is made to coincide 
with the Constant value of the potential in the region outside the spheres, 
then the integrals in (2.8) and (2.12) receive a nonzero contribution only 
from the part of the unit cell within the sphere. A trial function of 
the form, 


t l 

max. , ■ , 

<E> = S ? 1 
1=0 '' 


( 2 . 13 ) 


is used within the inscribed sphere, Rj^(r) is a radial function which is 
finite at r=0 and satisfies the radial differential equation 
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X i- (^2 + “i4) + v<r) 

2 dr . dr 2 ^ 

r r 


S 


Rj^(r) = 0 


(2.14) 


The functions lattice harmonics^ which are linear 


combinations of spherical harmonics of angular momentum 1. They are 


normalized, real, and mutually orthogonal. Finally, C are coefficients 


to be treated variationally . 


Then, Kohn-Rostoker have shown that the variational principle leads 
to a set of linear homogeneous equations 

1 

lax. ^ 

(2.15) 


Z 21Ah i-i- . 
l'=0 j* J ^1 


= 0 , 


y\ 1 j I’j' given explicitly below. The condition that the set of 

equations (2.15) have a nontrivial solution is the one that determines the 
value of E appropriate to the given k ^or value of k along, say, a line, 
for given namely 


det 


A 




I 4 I 


= 0 


(2.16) 


If in the expansion (2.13) of the trial function, 1^^^^^ is taken sufficiently 


large, the root of (2.16) approaches arbitrarily c'.o&e to the desired 
eigenvalue E = S(k) of (2.4). In practice, this convergence is extremely 


good, and it usually suffices to have = 2, 


The expansion of ^ ^ E^'O and r is given by 

^(r,r') = II ^ 

1 > j 1 j j 




(2.17) 
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where K = /\/e, spherical Bessel function, and Neumann function. 
yVn I'i' given in terms of B as: 

JsJ 

^ lj,l' j' ’ 


where we define 


r dR^(r) dj (iCr) 

*1 R.(r) dr j-(Kr) dr 
V 1 *^1 


Similar definition holds for 


i.rj' '[hj.i'r 


’ 1 

U' Qi<3i)J’ 


( 2 . 18 ) 


(2.19) 


( 2 . 20 ) 


where B^. are "Structure Constants", 

i j > i j 

The above expressions are appropriate for E^I^O, this zero of the 
energy scale is with respect to the constant value of the muff in- tin 
potential. The expressions for the case of E4^0 are given in Ref. 9. 

The same reference also contains the details about the structure constants. 


One great advantage of the Green's function method is that the 
structure part and the atomic part appear separately in the elements of 
the secular determinant. Apart from a scaling factor, the structure 
constants, constituting the structure part of the elements, are independent 
of the lattice constants. Therefore, although the structure constants 
are fairly ccmplicated functions, they may be computed once and for all 
for a given crystal structure. At the symmetry points for which the 
expansion (2.13) for 1 =2 admits only two or three independent terms, 

the actual evaluation of the (2x2) or (3x3) determinant of the 
'^Ij I'j' determination of the root of the secular equation (2.16) 
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may in fact be carried out with the use of a desk calculator, once 
structure constants are available. Alternatively, with availability 
of high speed digital computer, one can of course compute the required 
structure constants directly as a part of each run. 


In the application of the Green* s function method, therefore, 
one has to evaluate the structure constants for a sequence of values of 
E for given k, or for a sequence of values of k for given S. Next, the 
spherical Bessel functions and Neumann functions and their derivatives 

^ dRj(r) 

are obtained. The logarithmic derivatives, 


R^(r) dr 


1 


are obtained by numerical integration of the radial equation. The 

determinant of is then calculated, and the root of the secular 

1 j , 1 j 

equation (2.16) is found either graphically or by interpolation schemes. 



CHAPTER m 


ELASTIC SCATTERING OF ELECTRONS 
FROM SEMI-INFINITE CRYSTALS 


3.1 Introduction: 


A fundamental problem of the theory of electron-dif fraction is 
elastic scattering of electrons from the semi-infinite perfect crystals. 
This problem of elastic scattering from the semi- inf inite crystals arose 
from and is very closely related to what is at present the vast and 
complicated subject of LEED |Low-Energy Electron-Diff ractio^ problem. 


In 1928, a theory for LEED phencraenon was first proposed by 
14 

Bethe. His historical approach was an idealization of LEED problem 
as elastic scattering of the incident electrons from the periodic crystal 
potential. Then in principle it was a matter of matching the wave func- 
tions and their derivatives at the crystal surface. But this in actual 
practice involves lengthy computation before any reasonable result can 
be extracted out of it because of the slow convergence of the plaaae wave 
expansion for Bloch waves. The convergence difficulty was circumvented 
by Boudreaux and Heine^ who suggested the use of a pseudopotential 
instead of the crystal potential in Bethe' s matching formalism. They 
also pointed out the need to add evanscent waves to the superposition 
of Bloch waves used by Bethe. Similar type of formulation was also 
developed by Marcus and Jepsen.^^ All these techniques have certain 
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inherent difficulties like detailed calculation of the wave field inside 
the crystal and proper matching of the wave functions and their deriva- 
tives at the boundary. Various other approaches ’ " were proposed 

incorporating more realities for the LEED problem like surfaces etc. and 
using different approximations. But there also required heavy ccanputation 
and the computational difficulty continued. Recently, some elaborate 
and extensive computations ’ ~ were carried out on LEED, where the 

Important inelastic effects for LEED were accounted for, and led to very 
good agreonent with experiments. Still the basic computational aspect 
of the earliest and simplest idealization of LEED, as elastic scattering 
of electrons from the semi -inf inite periodic lattices, remained unchanged. 

The reasons, why one should pay attention towards the computational 
aspect and the detailed understanding of elastic scattering from the semi- 
infinite perfect crystals, are as follows: (1) Presently, the exact force 
law between an electron and a crystal is mostly unknown, and there is need 
to study elastic scattering in order to understand which of the predictions 
from such considerations of elastic scattering are sensitive to the various 

aspects of a force law describing the electron-crystal interaction. 

25 26 

Calculations in this spirit have been done before. * (il) Experimen- 
talists usually find it difficult to get into heavy computations to explain 
the observed intensities. They would prefer computationally simple and 
quick methods giving fairly good results to the complicated ones producing 
very accurate results, (ill) Understanding of elastic scattering in 
details is very helpful to have a clear-cut distinction of elastic effects 
and those due to Inelastic scattering -and surface effects. 
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Here the Green’s function method of energy band theory is properly 

formulated to be applicable in electron-diffraction theory* The KKR type 

3-7 

foimalisms for LEED problem are already there in the literature. But 

8 

we describe an approximate variational treatment of the problem which 

leads to satisfactory results without involving much conputation subject 

to the availability of standard programmes for band structure calculation 

5 

by KKR method. We consider the model like Kerre and Phariseau and 

others which assumes muffin- tin potential at the lattice sites of the 

semi- inf inite crystal and treats the problem variationally . General 

variational procedure for one-dimensional case is given by Shen and 
27 b 

Krieger, and Shen has also considered the case of the three-dimensional 
semi-infinite crystal. The virtue of the variational approach is that 
neither exact boundary between the crystal and the vacuum nor the detailed 
numerical value of the wave field inside the crystal are needed for the 
intensity calculation. All that is needed is a suitable trial wave 
function and the actual amount of ccmputation involved depends on the 
size of the determinants that one has to evaluate. However, there 
is no special computational advantage in a layer-by-layer variational 
calculation over the more usual ones by other methods. Moreover, the 
Green's function of such approach leads to fairly cumbersome structure 
constants . 

The above-mentioned difficulties of the variational approach are 
bypassed here by choosing Bloch wave as trial function. The wave vector 
of the Bloch wave is also treated variationally so that the approximation 
is more than a mere single Bloch wave approximation. Actually this amounts 
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to picking up the largest contributing Bloch wave from a combination such 
that its wave vector serves as a variational parameter. The choice of 
such an approximate trial function also finds sufficient justification 
in the final results obtained from it. It is shown that such trial 
function leads to the usual structural Green's function of KKR method 

9 

and then using lattice harmonics it is possible to show that the 
Intensity is a ratio of two very small size determinants whose elements 
can be formed out of those of the KKR secular determinant. Though the 
KKR type formalism is not new for this problem it never yielded to 
computation as simply and as quickly as the present one. However, the 


present approach will have the limitations of being applied only to 
the case of elastic scattering |*E-LEED^^^ intensities because of the 
constraint of using a static real potential. 


3*2 Variational Principle for Scattering Problems: 

The variational methods are very useful especially when perturba- 
tion methods become tedious and complex to use. Usually a single physical 
entity, in which one is interested, is given by a variational principle 
that leads to a stationary expression. Because of the stationary charac- 
teristic of the variational equation and its consequent insensitivity 
towards the errors in the trial function, it is often possible to obtain 
very good estimates for the physical entity under consideration with a 
relatively crude trial function. Obviously, this is a property of great 
significance in practical calculations. 
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Variational principle has been quite successful for the binding 
energy in atomic physics, scattering phase shifts in nuclear physics, 
and band theory in solid state physics. It ir also well known that the 
variational method can, in general, be applied fruitfully to the potential 
scattering problems. We describe, in this chapter, the use of variational 
principle for elastic scattering from (a) one -dimensional potential step, 
and (b) three-dimensional semi-infinite perfect crystal. In the following, 
we state the variational principle for reflection anplitude for one- and 
three-dimensional scattering problems. 


Consider the case of elastic scattering in one-dimension from a 
potential V(x) that is non-vanishing for The corresponding 
variational principle for the reflection coefficient R is 


= 0 , 


(3.1) 


and R is given by 
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R = (-^) 
o 


-j ik X -r ^ 

|Je ° V(x)'^(x) dxj 

J dx''p(x) V(x)£^(x) - J G(x,x' ) V(x? ) Y^x' )dx'^ 


wheTe'\j^(x) is the trial wave function, k^ is the wave number of the incident 

2 

plane wave, exp(ik^x); k^ = E, the total energy, and the Green's function 
G is given by: G(x,x') = -(i/2k^) exp (ik^ 1 x-x't ). 


, (3.2) 


For three-dimensional scattering problem, the variational principle 
for the reflection amplitude T which is a function of scattered wave vector 
k and Incident wave vector k. is 


0 , 


(3.3) 
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where 


ik, .r 


■ik .r 
-s 


(3.4) 


I? ( ’X — Ifr " — 'X ) ~\ 

I i/yA*(x)u(x)e dr )e ® (^^X^ r^^X^ <^X 

T(k ,k ) = — sssXsal — I == — (2 4) 

^ ^ J'y^*(r)U(r )j^'^(r)- ^ G(r,r ' )U(j;2' ''^^(x' ^“^X^ dr 

The ci'T vanishes for arbitrary variations of . Here G(r ,r ' ) is the 
Green's function satisfying (2.6), U(r) is the potential in atomic units. 


and V(r) is the trial wave function. 


3.3 Elastic Scattering From One- Dimensional Potential Step: 


We consider here, as an illustration, the problem of elastic 
scattering of electrons from one-dimensional potential step using 
variational procedure. In expression (3.2), it is possible to expand 


'f'(x) 


in terms of a complete set of functions as' 


(x) = ^ C. 0 (x), 


(3.5) 


and choose variations lly to get the condition determining R: 


i O i 


where 


=^0^(x)V(x) |^0^(x)- jG(x,x' )V(x' )0^(x' )dx'J dx, 


. ik X 

— J ° dx. 

For one-dimensional potential step 


(3.6) 


(3.7) 


0, X 0 

V^, X >0 


(3.8) 
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where is some constant value, and our choice of trial wave function 


Aj/' 


(x) = C(k) e 


ikx 


(3.9) 


with k in general ccxnplex; we have 

iV 


and 


D. 


11 


R 


(k + k )J ’ 

° 2 

iV (kk +k -V ) 
o o o o 

2kk (k+k ) 

o o - 


V k 
o 


(3,10) 


fk^ + 2kk ^+(k^-V )k -V kl 
to o o o o J 


If we want R to be extremum with respect to the variational parameter 
k, this extremum condition leads to: 


E = E(k) = (k^+V ) (3.11) 

o 

This is nothing but the energy wave vector relation for the case of 
one-dimensional potential step. 

% 

For E k = ijif, where (V^-E)^, which leads to the wave function 

-l< X 

e , and the reflected intensity is 

jRf^ = 1. (3.12) 


For E k = (E-V^)^, which is real, leads to the wave function 

xlcx 

e , and reflected intensity is 



(3.13) 


where, p* = k/k^, is called the index of refraction. Results (3.12) and 
(3.13) obtained here without using the usual matching procedure are the 
same as obtained by the procedure of matching at the boundary. 
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The important conclusion from this consideration is that the 

extremum condition for R leads to proper energy wave vector relation. 

For periodic potential the extremum condition would mean the energy 

27 

band structure relation. Actually, this is intimately related to the 
fact that energy band structure essentially governs the diffraction. 
Conversely, one can obtain the usual results variationally by choosing 
a wave function like (3.9) with k satisfying the known energy wave vector 
relationship. Thus R as a function of k is extremum for those values of 
k that satisfy the band structure relationship for the given energy when 
the potential is periodic . 


3.4 Scattering From Three-Dimensional Semi-Infinite Crystals : 

The simple model of a semi- inf inite perfect crystal is as follows. 
The lattice extends from - OO to +00 in x and y directions, and from 0 to 
+ 00in z direction. It is assumed that the crystal is a "simple crystal" 
having one atom per unit cell. Only elastic scattering is considered here. 

The reflected intensity in case of three-dimensional scattering is 


I = {Tl^/ien^ , 


where T is given by the relation (3.4) . 
For elastic scattering 


and here 




= k. 


2 

k = 


(3.14) 


(3.15) 


E. 
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E is the energy of the incident electrons. For diffraction from the 
perfect three-dimensional crystals, the potential is given by 

U(r) = 2 ^3*16) 

s 

where V(r-jr^) is the potential at r due to the ion at the lattice site 
r . In case of the semi- inf inite crystal, r lies in the half space 

o 

of 2^0. Then U(r) given by (3.16) is not periodic in translations 
along z-axis. 


Now, it is a matter of choosing an appropriate trial function '^(r) 
Note that it must have Bloch form along x and y axis due to periodicity in 
those directions. 


We assume the wave function 
following property 



within the crystal to have the 





(r+r ) 
s 


(3.17) 


for r and r in the region z^O. The crystal momentum is denoted by 

— — S r 

p. Thus, we are choosing a trial function which is exact for the x and 
y directions and which essentially contains the variational parameter p 
l^-component of pj. This would lead to the evaluation of the intensity 
as a function of p_ which has to be extremized to get the actual intensity 
Then 


T(k„,k,) = 





i<p-k.)-r f-7 


(k ,k.), 
—s ’’-i ’ 


(3.18) 


where 
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Here integration is over one unit cell and N is the total number of atoms 


in the crystal. Thus, one can evaluate T as a function of p and it is 


seen from the consideration of the one-dimensional problem that 


c>T 

"SpI 


= 0 


would give the p corresponding to energy E which may be used in T(p ) to 


evaluate the intensity, 


Thus, our evaluation of the intensity has all the virtues of the 

variational method. However, instead of taking a linear combination of 

Bloch waves, we have employed a single Bloch wave, the z-coraponent of 

whose wave vector is treated as a variational parameter. This is essentially 

the limitation of our calculation. But it is comparable to the limitations 

of the treatment of LEED by dynamical theory with a single- Bloch-wave 
29 

approximation, or the usual two or three beam band matching calculations, 

the s-wave multiple scattering theories, and the second order perturbation 

calculation. At tho s.ame time, the simplicity 'ind computational ease with 

our final expression is remarkable for numerical work. Much of the 

6 

simplicity of the computation will be lost if a linear combination of 
Bloch waves is chosen as the trial wave function. 


3.5 The Green's Function; 


Let us consider the Green's function in equation (3.19): 
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1 

N 



ip . ( r , - r ) 

Q(r4r «r* --r , ) e ^ ^ 

s “■ *^8 * 



S ip -ie** 

“£e»» ^ ^ 

tl ® 

SjS 

i£-Xs.. 

G(r ,r' +r „ ) e 
— ’— — s" 




where r^,, = r^, - r^. 

Reduction of the double sum to obtain the result in the third 
step involves the use of the actual form for G(r^,r^,) and the approKima- 
tion that the terms of the order of 1/N can be neglected. This can be 
seen by considering the one-dimensional case for simplicity [[see Appendi^ . 


It is easy to see that r „ is no longer confined to the half space 

8 

z^O. In fact, it is over the entire space and therefore is the usual 
KKR Green's function.^ Hence, the expression for_P"(k ,k ) becomes: 






ik,.r 

(r)V(r) e ^ "d^ 


'I 


-ik_ -r 


V(r) '^(r)dr ] 


] 


, ( 3 . 20 ) 


r , 

'f'* (r)V(r) |^^(r)-|(^(r,r')V(r' ^^(r' )dy dr 


with representing the usual KKR structural Green's function used 

in band structure calculation. This is not surprising because there 

exists an intimate relationship between the elastic scattering formalism 

4 

and the KKR method of band theory. In fact, the work of Watts establishes 
this connection and its use in the LEED theory. Here the variational 
version of this kind of approach is considered. 
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3 .6 Structure Factors and Bethe Condition : 

' It is very clear from the structure of relation (3.18) that the 
scattering intensity is proportional to 

^ I \ I 1 I \ / iCp-k ) *r . I 2 
N 


2 i(k.-£).rg 

^ 1 

"" N 

1 

s , 


1 ® ' 


which equals 


N 


2/3 


00 




i(k. -p )r 
i ^2 f 


1 e 


r =0 
s 

z 


00 


7 i(p -k )r , i 2 

' * z s s'* 

z z 
e 


(3.21) 


r , =0 
s ' 

z 


when 


II ®ll 


( 3 . 22 ) 


and zero otherwise. 


Here the symbols" || " and "z”, in the subscripts of the vectors 
kj,, kg and g j^a reciprocal lattice vecto^, denote the parallel and 
perpendicular j^i.e., z-componen^ components respectively to the crystal 
surface... of these vectors. This is because the summation over x and y 
components of r and r , yield delta functions which vanish unless the 

“*“S **~s 

condition (3.22) is satisfied. The condition (3,22) is nothing but the 

usual "Bethe Condition" that has emerged as it should from the structure 

factor parts of the variational calculation. This is comparable to the 

4 21 

work of Watts and Beeby and arises from the two dimensional periodicity 
of the lattice parallel to the surface. 

The two factors of (3.21) depend upon the particular crystal 
structure, and for cubic crystals one gets the following results: 
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(i) SC: 


r = La, 
s ’ 

z 


^ i(k^^-p^)al4 2 j 

1 

1=0 1 

0 

4 am %tk. -p )a. 
L 1 *^z 

z 


(ii) BCC; 


■s^ = hla. 




1=0 


-P.>T 

Z 


al 4 2 


4 sin %(k^ -p )a 


, (3.23) 


(iii) FCC: 


00 


s = %(l+m)a, 
z 


1 y i(k, 

CM 


f A > „ z z 


' 1 

1 r =0 
s 

z 

: 

16 sln^ ^(k. -p )a 
z ^ 


where 1 and m are zero or positive integers, and "a” is the lattice 
constant . 

3.7 Muffin-Tin Potential and Evaluation of Intensity ; 

We restrict our discussion to a potential, V(r), of the muffin-tin 
form. It is then possible to expand the trial function as 

1 

r-7 1 A-f 

^ - " ' ' " ^ (3.24) 


^£)= U c„<p> «1<P> 

1=0 j ■’ 


within the Inscribed sphere. 

We also expand, 

Ik^ .r 


^ *= 4p^ ^ i^ j^(kr) 


(3.25) 
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where and (0 ,0) are polar angles of and r, respectively. 

Now 



ik,-'; f t~ li- -r ^lf£ <r) "S 

‘x)V(r) e ‘ ar = ^ L e 


surface 


dS 


(3.26) 


Here use is made of the Green Theorem in vector analysis . 


We define D, .(k. ) by 
Ij -i 




and similar definition holds for D, .(k ). 

Ij -s 

Thus 


(3.27) 


ik^ .r 


(jr)V(i;) e ~ dr 


" S Di/i,) C*j(p) 

^ 5 3 


(3.28) 


Similarly 


I n ~ *— 

'^(r)V(r) e dr 


= 2 -^ D. , ., (k ) C , , (p) 

■*•^3 ° ^ 9 J ““ 

n 1 


(3,29) 


Thus the numerator of the relation (3.20) can be written as 

/ A D,.(k.) D,,.,(k ) C* (p) C,,.,(p) 

1 aTTV .. IJ -1 1 3 -s Ij r Vj' ^ 

^ 9 3 f ^ >3 

Similarly the denominator of (3.20) can be shown to be 

jp 


(3-30) 


J dr ( r ) V( r ) 

vol . or *“ jp 

m . t . sph . 


r ^ z,/^ ^^(r,r') ^ 

I C 

surf .of ^ i- 
m*t .sph . 


3 r‘ 4 


(r,r') -^fr— JdS 

at r*- 


= 2 F, ., C* (p) C 

1 4. 11 41 ^>J>^ >J ^>3 ~ ^ >3 

> 3 ♦ >3 


(3.31) 


If the zero of energy is adjusted to coincide with the constant 
value of V(r) in the region outside the inscribed spheres, the integrals 
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in (3.31) differ from zero only within the inscribed spheres. Within 
the inscribed spheres V(r) = V(r), i.e., a spherically symmetric 
potential . The integrals are over the volume and surface of the 
muffin- tin spheres. It has been already shown that ^(rjr') is the 

<5 ^ 

usual KKR structural Green's function. Therefore, the quantities 

are identical with ]_iji (2.18)|. Thus, denoting 

the combination (l,j) by n 



(3.32) 


30 j r 

Now, Kohn type transformations may be used to determine j^(k^,k^,p) . 

Thus, for the case D (k, )Tt D (k ) we define f (p) by: 

n ""i n "“S n — 


and 


fl^p) = S (3.33) 

n 

n 


to get 


H 

n,m 


F G*(p) 
n.m n *- 



^ V H f (p) f (p) 
n.m n A- m £. 

n,m 


(3.34) 


and 




) = 


f^(p) f2(p^ 


^ H f*(p) f (p) 

un^ n,m nP m?^ 


(3.35) 


where H can be expressed in terms of F , D (k. ), and D (k ) by using 

n^m * njui n X n b 


equations (3.33) in equation (3.34). 
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It follows that if ^ is to be obtained variationally , then the 
requirement that £7 is an extremum with respect to f^ is equivalent to 
the set of equations 


y [h 


2, ml ^n 


(3,36) 


which has non-trivial solutions if 


det H 


Thus , 


n,m " £7 '^l,n ^ 2,m. 


P - f det I 

~ " det 1 H 

U ' n,m 


(3.37) 


(3.38) 


where det j i® cofactor of Note that (3.38) gives 

uniquely, as the ratio of two determinants for a particular value of p, 
for which pjj j^omponent of £ parallel to the crystal surface^ is fixed 
by the Bethe condition leaving p^ ^ccmponent of p perpendicular to the 
crystal surface~| free. This result (3.38) is valid for the one-Bloch- 


wave trial function. 


As indicated earlier, p is to be treated as a variational parameter. 
There are two distinct cases: 

2 

(i) When E = k lies in a band, one can find p^ by setting 
det |f = 0. This is done by fixing E and varying p^ till the 

determinant passes through zero because as illustrated by the one- 
dimensional example, this is equivalent to band structure condition. 
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(ii) When E lies in the band gap, procedure (i) will not yield 
a zero to the determinant. Actually, p^ is to be taken as complex. 

This is equivalent to the case of E ^ for one-dlmcnsional example. 

It is essentially a method of finding E versus p relation when p is 
complex. 

3 .8 Discussion ! 

The formalism presented here clearly shows the extent to which 
the variational treatment can be carried out and the simplification 
that can be achieved with it. It also shows the intimate connection 
between the band theory, i.e., the KKR method of band structure calculation 
and the LEED intensity calculation. This fact is very helpful in view of 
the computational facilities available for the electronic structure 
calculations . 

That no physical aspect of the problem is lost because of the 

approximations involved in the formalism can be seen by a comparison with 
4 

the work of Watts which is based on multiple scattering theory. Thus 

the expression for intensity contains the cf^-function sums over reciprocal 

lattice vectors and the Bethe condition follows as a consequence of the 

periodicity of the system parallel to the crystal surface in both the 

approaches, that of Watts and ours. The KKR type matrix elements also 

appear in the Intensity expression of Watts. But our intensity expression 

1 2 

involving exact KKR * matrix elements is more easy for actual calculation 


compared to his result. 
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If we consider the computational aspect of the problem, the 
present formulation is not a complicated one. We have shown that it 
is possible to express the intensity uniquely as a ratio of two 
determinants expression (3.38)*|. The size of these determinants 

depends upon the number of terms in the trial function |^given by (3.24^ . 
For low energy electrons, this number is certainly small and so is the 
size of the determinants. And if the symmetry of the problem is 
exploited fully, as has been done here by using lattice harmonics, the 
size of the determinants reduces further. Since the elements of the 
determinants in (3.38) are expressed in terms of exact KKR matrix 
elements, the numerical computation of the intensity would not pose 
any problem once the KKR band structure programme is ready. The 
advantage of this variational calculation is that it avoids matching at 
the boundary and the detailed calculation of wave field inside the crystal. 

Of course, the limitations of the approach come from the use of 
an approximate trial function pi single-Bloch-wave approximatioi^ , and 
a real static potential pi*®* 

muffin-tin form. This treatment is employed in the next chapter to 
compute the (0,0) intensity of normally incident electron beam, where 
we will compare our results with some other theoretical and experimental 


exclusion of inelastic effectsT ox the 


data and discuss the results. 



CHAPTER IV 


LEED SPECTRUM AND BAND STRUCTURE OF xlUU/.INIUM 


The formalism developed in the last chapter is applied here to 
compute the (0,0) intensity from (100) surface of Aluminium crystal. 

In the following section, we formulate the case of normal incidence 
and back scattering for which immense simplification occurs. Next, we 
describe our results for the high-energy band structure of Aluminium, 
and at the end the results are discussed. 

4.1 Formulation : 

For the case of normal incidence, we take the incident beam of 
electrons to be perpendicular to the crystal surface, and treat this 
as the z-direction. Then 

k^ = (0, 0, k), (4.1) 

0 = 0 . =0 and k , = 0 . 

1 i|| 

Now, we proceed to find out back-scattering-intensity; 

Here 

k =0, (4.2) 

®ll 


and hence 
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Thus in this case we have 

Dy (kj) = %(0,0) 

■= Oij(k^) , (4.3) 

p s (0, 0, p), 
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and the Kohn transformations described by (3.33) do not apply. 
However, in this case of = Dj^^(kg), one can define 


8i = S 


ind 


Then 


g = C(p) forni^l. 
®n n 




8 1 8 


S 8„ 8, 

n,ra 


m 


It follows that 


whence 


r<p) 


1 

r 

<^n,l t^m.l 


hi.il 


det f G j 


n,m| 

r- • 


n 


= 0, 


where det 


hi, 4 


is the minor of G 


1,1' 


Now to evaluate the elements G^^^ we note that 




(4.4) 


(4,5) 


(4.6J) 


(4.7) 


(p) = 


^ 1 . 1 1 • I s'f - Sii'i, (4.8) 
1,1 lj,l J ^Ij^l'j' 


lj,l' j 


and using (4.4) one gets 
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00,00 


00,1'j' 


lj,00 




00.00 
^ ^^00^ 4 


^00,00 i' 


’C? 


r',1’ ’* 


F „ D 
00,00 li 

2 


F + F 

^ ^ A 114* 


^00 

li,00 

^00 

D. .D 


11,111' 


- F, 




lj,l'j- 00,00 , ,2 00,1'j' D 

00 ^ 


- F,, 


D^. •■ I 1,1' ^ 0 




j,j' ^ 0 I . 


Now, substituting from (2.18) and (4.3) in (4.9), and after doing 
algebra, we get 


f detfm'l 

Plj.l'j' 


det 


r<p) 

where det | m' | is the minor of After putting 


and 




.( 0 , 0 ), 

the elements G' , are given by: 

i J , i J 

Aoo,oo ’ 


00,00 


‘^00,1'j' 


""JAoO,!' j'^00 ■ AoO,00 ^1' j'J , 


(4.9) 

some 


(4.10) 


(4.11) 
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••• (4,12) 


4.2 Band Structure Calculations ; 

It is already mentioned, that the knowledge of the real and complex 
band structure of a crystal is required before taking up the intensity 
calculations. The KKR method of band structure calculation is adopted 
to compute the electronic structure of aluminium. The choice is aluminium, 
because lot of experimental information is available for it. 

For the energy E in the band gap, the electron wave vector p is 
taken as complex. Then one specifies E within the band gap and solves 
for the real and imaginary parts of £ as function of E. For complex p, 
the Bethe condition stays, but structure facco.s end structure constants 
get changed. The expressions for structure constants obtained by Kohn- 
Rostoker^ are no longer useful for ccmplex £ case. In that case, one has 
to evaluate the structure constants ab initio. Kcrre-Phariseau have 
derived ejiplicit expressions for structure constants, when p becomes 
complex. 

There are a variety of problems, like LEED intensity calculations, 
which need to probe the electronic states of a solid well above the Fermi 
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FIGURE 1: ENERGY BAND STRUCTURE OF Al 
ALONG A -AXIS 
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Figure 1 shows 


level ^i.e., high-energy band structure of a solid 

our results for high-energy band structure of aluminium along the 

symmetry direction, A. Only few other calcu^ ^.tions of the high energy 

band structure of aluminium have been performed by Connolly ^self- 

T 32 r 

consistent APW calculation!, Hof fste in-Boudreaux £ OHW based pseudo- 
potential calculation^ , and recently by Jepsen et al7 [ftcm eigenvalues 
of transfer matrix^ 


33 

Self-^-consistent potential for aluminium obtained by Snow is 
employed. The programme was tested for and found to give the well-known 
rapid convergence of KKR method. We also used Mathieu potential and 

9 

found that the results of Hara-Segall are reproduced. The programme has 
been found satisfactory for the range of energy involved in LEED intensity 
calculations. The results in both low and high energy region are in fair 
agreement with the previously published data. In low energy region the 
results emphasize the free-electron like nature of aluminium. 


The energy bands also have been compared with the calculations by 

7 3 3 

Connolly^ ^ and Jepsen et al.' with the same potential for aluminium. 
Careful comparison shows that but for a few discrepancies in limited 
regions the results are very close. It is important to point out that 
any residual inaccuracy is unlikely to giffect the variational calculation 
of the LEED intensity. 


4.3 Results and Discussion; 


Ultimate justification of any variational calculation lies in the 
quality of the final results. In this spirit we have computed the (0,0) 



ENERGY (eV) ^ 

FIGURE 2 : COMPARISON BETWEEN MEASURED AND CAlCUtATED BEAM INTENSITIES 
VERSUS ENERI^Y FOR THE (0,0) BEAM FROM THE (100) SURFACE OF At 
NORMAL INCIDENCE. THE ZEROS FOR CURVES a,b,c MARKED WITH 


%xpt/ ^Th>^> RESPECTIVELY, 
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ENERGY (RV) 


FIGURE 3: COMPARISON BETWEEN MEASURED -AND CALCULATED 
BEAM INTENSITIES VERSUS ENERGY FOR THE (0,0)BEAM 
, FROM THE (100) SURFACE OF Al , NORMAL INCIDENCE . 
THE ZEROS FOR CUVES a,b,c MARKED WITH 
%xpt » Oib * ° ^ RESPECTIVELY. 
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intensity from (100) surface of A1 for normal incidence, as a test case. 

The Figs. 2 and 3 compare our results with some other experimental and 

theoretical results. In both the figures, the curves "a" are experi- 
,34,35 

mental curves, the curves "b" are obtained from band -matching 

29 3 6 

calculations, ’ and the curves "c’* are the computed results from our 
variational formulation. We have included, in these calculations, 
only the first few 1 values in view of the rapid convergence in 1 that 
one gets for the KKR method. The accompanying high-energy band struc- 
ture calculation of A1 by KKR method justifies this convergence for 
the energy range involved in the present calculation. The contribution 
estimated for higher 1 values is very small and may be ignored for the 
present purpose . 

While comparing the above curves one must keep in mind that our 
model is a gross simplification of the experimental situation of LKED 
and that the reproduction of the fine structure, exact peak heights 
and widths ‘^of the measured LEED intensity curves has proved to be an 
extremely difficult task even for the most elaborate computation. 
Therefore, ignoring such "fine" effects, there is a great similarity 
between the experimental curves and the calculated curves. CXir results 
in Fig. 2 are better than those from band matching calculations with a 
"Single-Bloch-Wave-Approximation '* while in Fig. 3 they are comparable 
with Hoffstein-Boodreaux calculations. With the possibility of 
numerous other physical realities being brought into the calculation, 
such an agreement is remarkable and provides ample justification for 
the chosen trial wave function. 



CHAPTER V 


ELECTRONIC STRUCTURES OF BERYLLIUM WD TITANIUM 

5.1 Introduction : 

Electronic structures of Beryllium and Titanium, the two typical 
hexagonal close-packed metals, are of considerable theoretical 

and experimental interest. In literature, there exist theoretical 
calculations ’ of their band structures using methods other than the 
Green's function method. We report, in this chapter, our theoretical 
results^^ on electronic band structure of Beryllium and Titanium obtained 
by employing the Green's function method. 

As already mentioned earlier, the Green's function method deserves 

wide application. But the method has found application mostly to cubic^° 

11 

metals and limited application to other structures. There are two 
reasons for this, namely, the computational difficulties encountered 
while computing the band structure of complex crystals and the inadequacy 
of the muffin-tin approximation in case of structures which are not cubic. 
However, the case of the hep metals is different, because for the hep 
metals the muff in- tin approximation is almost as good as for the cubic 
metals. It is also possible to simplify the computation to a great 

13 

extent by the careful application of group- theoretical factorization 
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of the secular determinant* Then it is interesting to apply the Green* s 
function method to compute the band structure of the hep metals and fully 
exploit its virtue of fast convergence of eigenvalues. This has been a 
major consideration of the present study which also provides useful 
information for comparison with the results hy other methods. 

Beryllium, the lightest of all the hep metals, has a simple atomic 
2 2 

configuration (Is) , (2s) . Its band structure is very well known and 

the first theoretical band structure calculation^*^ was done as early as 

in 1940. For the present study, its serves as a test case. On the other 

38 

hand, the band structure of Titanium has been studied rather recently. 

It is a typical transition metal with more accepted atomic configuration 
2 2 

(3d) , (As) . In view of the limited experimental information and a few 
theoretical calculations, its band structure needs to be investigated 
further. The Green's function method is a valuable tool to carry out 
such investigations. The results are expected to exhibit the general 
features of the hep transition metals. 

5,2 Brief Account of the Method ! 

Green's function method for "conplex" cryTtils was developed by 

1 o 2 

Segall, in a manner parallel to that of the Kohn-Rostoker method for 

"simple" crystals. In this section, we give just an outline of the 
1? 13 

method and describe the group- theoretical factorization of the 

41 

secular determinant employing hexagonal harmonics for symmetry point* 

and symmetry directions of the Brlllouln zone. The details of the 

12 

method have been discussed at length in the literature. Here we omit 
any such detailed discussion. 
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The complex crystal is defined by the lattice translation vectors, 

r^, given by r^ = n^a^ + n 2 a 2 + n 3 a 3 ^n^ = 0,+l,+2, where a^, 

£-2 ^3 basis vectors of the lattice; and the set of vectors 

~1’ ~2’ '*■’ — which locate the centres of ^ atoms in the unit cell. 

Variational principle as in (2.11) is set up and a procedure very similar 

2 

to the one by Kohn-Rostoker for the simple crystals is followed to solve 
the problem of band structure calculation for the complex crystals. The 
important difference between the two procedures ’ is that the quantity 
i/\for the complex crystals []]see (2.12^ contains terms with two types 
of products of surface integrals unlike the case of simple crystals. The 
first type of products arise from integrals over surfaces centred about 
the same ion |^terms to be called "diagonal" in the centres^, and the 
second type of products are due to integrals over surfaces centred about 
different ions j^erms to be called "off-diagonal" in the centres^. This 
is a result of having more than one atom per unit cell in the complex 
crystals. The terms off-diagonal in the centres are slightly more 
troublesome to evaluate because the natural coordinate systems for the 
position vectors r .and r' in the muffin-tin spheres centred about two 
different ions have their origin at different points. But after 
redefining the position vectors and the structural Green's function 
the off-diagonal terms are readily evaluated. Thus, the quantity 
becomes the sum of all diagonal and off-diagonal terms. Finally, 
simplifying the whole thing one arrives at the band structure condition 




( 5 , 1 ) 
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where 


jtLlm.V 


m 


K 11 1 + 

Im, 1 m 


Q,f\^ 1 

Q^( j^) J r=rj, 


(5.2) 


Qi<li)-| 

* The expression, -— ” j-j j evaluated at r = , the muffin-tin 

1 1 - 1^2 
radius of the sphere in the unit cell. Structure constants 

(II \|it \ 

A, ,, ,, which are diagonal in the centres, are identical to the 
Im, I'm 

2 

structure constants A , , in Kohn-Rostoker treatment for a simple 

Xm X in 

crystal with the same translational symmetry. Whereas the structure 

constants off-diagonal in the centres, different. 

12 

They are discussed in details by Segall, 


The secular equation (5.1), while being valid for all k (^crystal 
momentum vect<^ vectors including those with special symmetry, can be 
considerably simplified at symmetry points and along symmetry axes as 
follows : 


Consider simple lattices with only one atom per unit cell for 
which the trial wave function within the muffin-tin spheres is given by 

■(r) =. Z C^.RiCr) 

I 




(5.3) 


where R^(r) radial solution of the Schrodinger equation, « 
the spherical harmonics and are variational coefficients. One 
use symmetry-adapted functions 


can 



z 

m 


v*^ Y, (e, 0 ), 

1 ,m Im 


(5.4) 


by determining the coefficients of the unitary matrix V* group- 

theoretically. Using these functions, obtains the 
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secular determinant^ which can be related to A,, ,, , obtained 

J lm,l'm 

by using ^^^(9,0) through the relation 





£jk 



A 


Im, 1 'm' 


j 

' ,ni 




whence (5,5) 

A A A A 

A = v,A.vt . 


A 

As is well known, A will be in the factorized form ^pr block diagonal 
form^. Extending the same idea to the case of hep structures with two 
atoms per unit cell, one can easily factorize the determinant. Now the 
trial wave function can be written as 




RjCr) Yi„(§. 8 )], 


(5.6) 


where quantities with bars at the top refer to the second atom in the 
unit cell as the origin. Again the symmetry- adapted hexagonal harmonics 
may be used to affect factorization. It is Important to point out that 
if hexagonal harmonics are used directly in the expansion of , one 

gets involved with two-centred integrals which are very conveniently 

A 

avoided by using the corresponding transformation matrix V. Thus, the 
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factorization using V, which is trivial in the case of one atom per unit 
cell is very helpful for computation involving more than one atom per 
unit cell. 


5.3 Description of the Calculation : 

The crystal potentials used for determining the energy eigenvalues 

*^7 38 

are derived by Terrell ' for Beryllium and Hygh-Welch for Titanium. 
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I 

FIGURE 4: THE B.Z. FOR THE HOP CRYSTAL | 

WITH THE 1/24 ZONE AND CRYSTAL | 

SYMMETRY DIRECTIONS. i 
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TABLE I: Points of 1/24 of the Brillouin zone x<ihere energy eigen- 
values are computed 1 Wave vector k = 

% i 1 24 2 


'Y ^ 2 ^ 1 i A 

+ bo , where = (-=;—) k , b„ = (-;^k +k ) and 

d| 1 a X 2 a X y 

1 A 

^3 ~ c basis vectors of the reciprocal lattice with 

A /\ A. 

a and c as the lattice constants, k s k and k are unit 

X y z 

vectors along x, y and z directions, respectively^. 

Label 


Weightage Factor 
(No. of like vectors) 


Wave Vector 

r 



1 



(0,0,0) 

s 



6 



(3,0,0) 

Z 



6 



(6,0,0) 

Z 



6 



(9,0,0) 

M 



3 



(12,0,0) 

T 



6 



(2,2,0) 

- 



12 



(5,2,0) 

- 



12 



(8,2,0) 

T' 



6 



(11,2,0) 

T 



6 



(4,4,0) 

- 



12 



(7,4,0) 

T' 



6 



(10,4,0) 

T 



6 



(6,6,0) 

T' 



6 



(9,6,0) 

K 



2 



(8,8,0) 

A 



2 



(0,0,6) 




12 



(3,0,6) 
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TABLE 1 ( . .Contd.) 


Label 

Weigh tags Factor 
(No. of like vectors) 

Wave Vector 

(^,9,ry 


12 

(6,0,6) 

- 

12 

(9,0,6) 

U 

6 

(12,0,6) 


12 

(2,2,6) 

- 

24 

(5,2,6) 

- 

24 

(8,2,6) 

- 

12 

(11,2,6) 

- 

12 

(4,4,6) 

- 

24 

(7,4,6) 

- 

12 

(10,4,6) 

- 

12 

(6,6,6) 

- 

12 

(9,6,6) 

P 

4 

(8,8,6) 

A 

1 

(0,0,12) 

R 

6 

13,0,12) 

R 

6 

(6,0.12) 

R 

6 

(9,0,12) 

L 

3 

(12,0,12) 

S 

6 

(2,2,12) 

- 

12 

(5,2,12) 

- 

12 

(8,2,12) 

S' 

6 

(11,2,12) 
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TABLE I ( . .Contd . ) 


Label 

Weightage Factor 
(No. of like vectors) 

Wave Vector 

<■(,? ,r) 

S 

6 

(4,4,12) 

- 

12 

(7,4,12) 

S' 

6 

(10,4,12) 

s 

6 

(6,6,12) 

S' 

6 

(9,6,12) 

H 

2 

(8,8,12) 
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These potentials were calculated by superposing the self-consistent 
Hartree-Fock atomic potentials and introducing the exchange by means 
of Slator s free-electron exchange approximation. Apart from the fact 
that these potentials are satisfactory and carefully constructed, our 
reasons for using them Is also to ease the comparison of the resultant 
electronic structure. 

Figure 4 shows the Brillouin zone for the hep crystal along 
with the (1/24) of the zone used in the calculation and the various 
points and axes of symmetry. The hep metals contain two atoms in the 

12 9 

unit cell so that one has to compute off-diagonal structure constants * 

which are, in general, complex. Energy eigenvalues are computed at 

45 points in (1/24) of the Brillouin zone. The coordinates of these 

points along with their welghtage factors are given in Table 1. Using 

41 

the symmetry-adapted hexagonal-harmonics up to 1=2 we obtained reasonable 
convergence of eigenvalues. 

The Fermi energies, E^, are obtained by using the graphically 

interpolated energy eigenvalues which are equivalent to 1152 points in 

the Brillouin zone, as follows: There are four electrons per unit cell 

for Beryllium and eight electrons per unit cell for Titanium and each 

band can accommodate two electrons per cell. Thus, for a discrete 

mesh of 1152 points, there will be 2304 energies below Fermi energy for 

Beryllium and 4608 energies below Fermi energy for Titanium. Arranging 

the 2304 energies for Beryllium and 4608 energies for Titanium in 

ttl 

ascending order and counting up to the 2304 energy for Beryllium and 
4608*"^ energy for Titanium yields the corresponding Fermi energies. 




VECTOR, A.U. 

5 IN (b) DOUBLY DEGENERATE. 
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The counting scheme, mentioned above, is also used to calculate 
the density of states histograms. This is carried out simply by counting 
the number of the ordered energy eigenvalues that fall within each of 
a fixed set of energy intervals which cover the calculated energy range. 

A number of histograms are constructed by varying the energy interval 
or bar width, up to the stability situation when no appreciable 

change takes place in the histogram due to a small alteration in^E. 

For Titanium, the following study is also carried out at a 
preliminary level. First, the band structure of Titanium is calculated 
for a value of the exchange parameter, , lying in the range of 2/3 
to 3/4. This is interesting because it was concluded from an earlier 
work^® that 2/34(<’^3/4 is the possible range of the exchange parameter 
that could produce results close to reality. Next, the effect of varying 
the atomic configuration on the band structure of Titanium is studied. 

This is interesting from the point of view of proper interpretation of 
experimental information available regarding the band structure and it is 
closely related to the well known Inherent ambiguity for transition metals. 
In both the cases, the crystal potential was constructed by superposing 
the self-consistent Hartree-Fock potentials, and the exchange was treated 
in the Slater's free-electron exchange approximation. 


5 . 4 Results and Discussion : 


A. Band Structure : 

The electronic energy bands for Beryllium al ong the sy mmetr y 
directions are shown in Fig. 5. We have shown the 


je ir s ^ . thjeej^ihds fj^r 

CENTH/\L 


Acc. No. 




3082-7 
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Occupied Bandwidth, 

L®F"F^3» Rydberg. ^ 


Occupied Bandwidth, 

Rydberg. 


Theoretical: 


Theoretical ; 


Herring-Hill 0.865 

^ 42 

Jacques 

Single OPW 0.956 

OPW with variational 

Parameter 0,919 


Hygh-Welch' 


Se If - c ons is tent 
APW 

Present 


0.211 


0.153 


0.240 


Cornwell 


Loucks-Cutler 


Terrell' 


Present 


0.840 


0.780 


0.901 


Experimental : 

i 

Skinner et al . 


0,84 + 0.01 


0.868 


0.441 
(6 eV) 


Experimental : 


Skinner 


1.01 + 0.07 



54 


Beryllium. Since there are four atoms per unit cell and each band is 
doubly degenerate, on the average there will be two filled bands for 
Beryllium. The band structure of Titanium along the symmetry directions 
are shown in Fig. 6. There are eight atoms per unit cell in this case 
and on the average there will be four filled bands out of the first 
five bvands depicted in Fig. 6. 

37 

A comparison between our results and those of Terrell by the 

Augmented Plane Wave J^APwJ method for the case of the Beryllium shows a 

great deal of similarity. Since the electronic structure of Beryllium 

is well understood, this shows that the present application of the 

Green's function method has proved to be sound in this test case. Now 

38 r*^ 

a comparison of the results of Fig. 6 with those of Hygh-Welch APW 
calculation^ shows that they are very similar. 


B. Fermi Energy, Ej,: 

Fermi energies are obtained in both the cases by the counting 
scheme described in the previous section. We obtain Eg, = 0.83 Rydberg 
for Beryllium and = 0.335 Rydberg for Titanium. Table II compares 
our results with other theoretical^^ '^nd experimental 
results. While our results are very close to those obtained by the APW 
method in both the cases, the experimental results are significantly 
different. As usual, a consideration of the difficulties inherent in 
the extrapolation of the low-energy tail of the x-ray spectra leads to 
the conclusion that the departure of experimental values from the 
theoretical values is primarily due to this extrapolation. The limited 

derivable from the present calculation indicates 


Fermi surface information 
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no Significant departure from the corresponding APW work which had 
been rather exhaustive In this respect. 

C. Density of States, D(E) ; 

The density of states histograms are obtained by following the 

procedure outlined in section 5.3. The resulting stable histograms 

for Beryllium are obtained for Ae =0.1 Rydberg and is shown in Fig. 7. 

37 

It resembles very closely with that obtained by Terrell . The basic 

feature of the density of states curves remain the same in the work of 
40 44 37 

Herring-HlH, Loucks-Cutler, Terrell, and the present calculation 

It is also interesting to point out that the density of states 

curve for Beryllium shown in Fig. 7 is consistent with the observation 

47 

of Johnston and Toraboulian, that the x-ray absorption edge coincides 

with the high-energy limit of the emission line in Beryllium and that 

the absorption data have a peak at the low-energy end. However, an 

extension of the density of states calculation to a higher energy is 

necessary for further comparison with the x-ray data of Johnston- 
47 

Tombou llan. 

The deep minimum in the density of states curve for divalent 

elements like Beryllium is not the characteristic of the hexagonal 

structure. But this deep minimum has special relevance in the case of 

Beryllium because its metallic properties depend upon the overlapping 

2 

between s and p bands. The band up to Fermi energy represents (2s) 
electrons of the Beryllium atom and the density of states of course 
goes to zero in the gaps between energy bands. It is clear from all 


39 



ELECTRONS/ATOM-RYDBERG, D(E) 
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TABLE III, Heat Capacities 


Be 

Electronic Specific 
Heat Coefficient,'/ , 
in units of 

Cal.mole^’ deg. 

Electronic Specific 
Heat Coef f icient , y , 
in units of 

-1 -2 

Cal. mole^ deg. 

Theoretical : 


Theoretical : 


37 

Terrell 

0.3x10"^ to 0.5x10"^ 

Hygh-Welch^® 


Present 

0.41x10“^ 

APW 

11. 8 1x10“'^ 



Present 

7.263x10"'^ 



Hygh- Welch®® 




Self-consis- 
; tent APW 

9.18x10“^ 

Experimental 

: 

Experimental ! 


Hill-Smith^ 0.54xl0“^ 

Daunt 

8x10"^ 

49 

s -4 * 

(O.44+0.04)xl0 1 

51 

Gschneider 

7.223x10'^ 




this that the deep minimum shows in a graphic way how the overlapping 
of the (2s) and the (2p) band is not really very great. 
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The typical density of states curve for hep metals is exhibited 

by Titanium as shown in Fig, 8, The histogram agrees well with the one 

38 

obtained by the APW method. There are two peaks due to the relatively 
narrow d bands with the Fermi energy on the high-energy side of a 
relatively sharp minimum. 

Theoretical estimate for low temperature electronic specific 

heat coefficient,^ , is calculated by determining the density of states 

at Fermi energy D(i^). Our results for both Beryllium and Titanium 

37 38 A8-51 

are compared with other theoretical ’ and experimental results 

in Table III. 

The estimates of '^f rom the present calculation may be seen to 
agree fairly well with other theoretical values in Table III. For 
Titanium the value of ^ is much higher than the free electron value of 
1.15 X 10~ cal .mole"! deg.” The d bands are responsible for the 
Increase in the specific heat value from its fror electron value. It 
may be expected that since the electron-phonon interaction has been 
neglected in such estimates, there will be some discrepancy with the 
experimental values. This is borne out by the figures presented in 
Table III. Titanium is seen to have much less effect of electron-phonon 
interaction than Beryllium which is very light ccmipared to it. 

D. Atomic Configuration and Exchange Earameter in Titanium : 

The reasons for studying the electronic structure of Titanium with 




WAVE VECTOR, A. U. 

FIGURE 9: BAND STRUCTURE OF TITANIUM FOR DIFFERENT ELECTRONIC CONFIGURATIONS AND DIFFERENT VALUES 
♦ OF EXCHANGE PARAMETER 
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different atomic configurations and different values of exchange 
parameter are outlined in the previous section. We have studied 
these aspects at the preliminary level which is perhaps sufficient 
to show the qualitative features of the changes that may be expected 
under these conditions. 

The results presented in Fig. 9(a) and Fig. (9b) correspond 

to the band structure of titanium along FT’K for the atomic 

3 1 2 2 

configurations (3d) , (As) and (3d) , (As) ; respectively. Exchange 

parameter, tC » is taken to be unity in both the cases and the crystal 

potentials used were obtained by superposing the self-consistent 

Hartree-Fock atomic potentials with exchange treated in the Slater's 

f ree-electron approximation. This kind of variation in the atomic 

configuration, in general, results in a narrowing of the (3d) band 

which is evident from the figure. It is admitted that such computations 

52 53 

cannot remove the Inherent undertalnty ’ of atomic configuration 
for transition metals. Perhaps the results of self-consistent energy 
band calculation coupled with detailed and unambiguous experimental 
information can lead to a satisfactory resolution of the question. 

Figure 9(b) and Fig. 9(c) represent band structure of Titanium 

2 2 

along pT K for the atomic configuration (3d) , (As) with exchange 
parameter .C = 1 and / = 0.7, respectively. The second value of|£ = 0.7 
lies in the interesting range of values anticipated from the APW 
calculation.^® In Table IV, we show the approximate estimates for 

occupied band widths “ £*1 1 ^ band widths ^ 

^ ,38 , . _ ,46,5A555 

them with other theoretical and experimental 


compare 


estimates . 
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This preliminary study indicates that / =0.7 leads to better agreement 
with experiment than/= 1. However, this is no definite conclusion 
in view of the limited experimental data on Titaniimi but it surely 
supports the previous conjecture regarding the interesting range of / . 
The exact value of will remain an open question for Titanium as 
long as detailed experimental information is not available. 



CHAPTER VI 


CONCLUSION 


Some interesting conclusions, drawn from the present study of 
"Electron-Diffraction" and "Band Structure" problems by the Green's 
function method, are as follows. 

The variational treatment described in chapter III for the 
calculation of elastic scattering intensities from a perfect semi- 
infinite crystal avoids the calculation of cumbersane structure 
constants by assuming a Bloch wave trial function. The layer-by-layer 
treatment of the problem appears unnecessary for the case of a perfect 
semi- inf inite crystal. Be the condition follows, as it should, from 
the structure factor part of the variational calculation. From the 
standpoint of computation, it is importaiit cj exploit fully the 
symmetry of the problem. This has been done by using lattice harmonics 
in the given variational formulation. It has been shown that ultimately 
the intensity may be found by evaluating the ratio of two very small 
size determinants which immensely simplifies the computation. And, 
it has been concluded that such variational calculation of intensity 
can be performed as a byeproduct of the band structure calculation 
and this fact is very helpful in view of the detailed computational 
facilities existing for the latter. 
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The comparison, in Chapter IV, of the computed intensity 
profiles with the experimentally observed spectra implies that, in 
practice, one can use the approximate variational treatment given 
in Chapter III and evaluate elastic LEED intensity curves which can 
agree within a few per cent with the experimental results. Of course, 
there are other elaborate computational schemes which lead to very 
good agreement with experiments. But they involve a lot of labour 
and ccanputer time. It all depends upon whether one ^e.g., an 
experimental is can devote much of time and effort for such computa- 

tions and whether all the detailed fine structures of LEED intensity 
spectrum are of interest in any particular case. The method of 
Chapter III is useful when the efforts that can be put on computation 
are rather limited. It may be a choice of quick and inexpensive 
calculation leading to approximate results over a time consuming and 
expensive calculation producing rather accurate results. This is the 
conclusion that follows from comparison of our results with those of 
LEED intensity measurements . 

The model considered in Chapter III for intensity calculations 
is a gross simplification of the experimental situation of LEED and 
to include more realities one has to consider the important effects 
like inelastic effects, surface effects, etc. But this would mean 
putting in tremendous efforts on computation. 

The results of electronic structure calculation for Beryllium 
and Titanium by the Green's function method provide an independent 
first principle check on the previous theoretical calculations, and 
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such efforts are in accordance with the best traditions in the Band 
Structure Theory.' This was not required in case of Beryllium for 
which quite exhaustive information is available. But the same 
cannot be said about Titanium. In fact, it is a candidate for still 
further theoretical and experimental investigations. Perhaps several 
others like Zirconium etc. in the same period as Titanium should also 
be studied to resolve the questions regarding atomic configurations 
and exchange parameter values • 

Out results of band structure calculation are quite close to 
APW results. This may stimulate interest for further theoretical 
and experimental investigations in case of Titanium. The Green's 
function method could be a handy tool for such further studies. But 
no doubt that more detailed and reliable experimental information 
will be needed before unambiguous theoretical answers can be provided 
for many Important questions. This is exactly why our inferences 
about for Titanium must remain at the level of a conjecture. 

Finally, we conclude that it is possible to obtain satisfactory 
numerical results for both the electron scattering problem and the 
electronic structure problem of the hexagonal close-packed metals by 
employing properly formulated and carefully simplified form of the 
Green's function method of Korringa-Kohn-^ostoker. 
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APPENDIX 


Assuming the lattice constant to be unity, one can define 


f ( s 


-£* ) = ^ G(r-r* "ts-s* ) e ® ^ 


Then, 


N N 


N N-s’ 


SI S£‘s-s') ' S S ,f( 

s=0 S“0 s -0 l=-s' 


1) , 1 = s-s' 


N N-s' N -1 

> n s * tL n £(i) 

s'=0 1=0 s‘=0 l=-s' 

N N-1 -N N 

= n * s n £(') 

1=0 s'=0 1=^ s'=-l 

N -i 

= 21^ (N-l+i)f(l) + 2^ (N+1+1) f(l) 

1=0 1=-N 

+N +N +N 

= E Nf(i) . H f(i) - E lu f(i). 

1=.N 1=-N 1=-N 


The second term is obviously of the order of 1/N and can be 
neglected compared to the first one. The third term is also negligible 
because 

ik^fx-x'l 

G(x-x') = (-i/2k^) e , 

ik ^l|-ipl 

so that for large 1, fd) iSi (-l/2k ) e 


Therefore, 
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and 


|I 111 £(i) = (-i/k,) S 


f (1) 


which is of the order of i/N of the first term. Thus the desired 
result follows. In case of three-dimension, x and y direction 
lattice sums are trivial and exact and the sum irv z-direction needs 
the above-mentioned approximation. 



